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TRANSPORT ANALYS IS OF MEASURED NEUTRON LEAKAGE SPECTRA FROM 
SPHERES AS TESTS OF EVALUATED HIGH-ENERGY CROSS SECTIONS 
by Donald Bogart, Donald F. Shook, and Daniel Fieno 

Lewis Research Center ' 

SUMMARY 

Integral tests of evaluated ENDF/B high-energy cross sections have been made by 
comparing measured and calculated neutron leakage flux spectra from spheres of vari- 
ous materials. An americium-beryllium (Am-Be (a, n)) source was used to provide 
fast neutrons at the center of the test spheres of beryllium (Be), polyethylene (CE^), 
lead (Pb), niobium (Nb), molybdenum (Mo), tantalum (Ta), and tungsten (W). The ab- 
solute leakage flux spectra were measured in the energy range 0. 5 to 12 MeV using a 
calibrated NE213 liquid scintillator neutron spectrometer. Absolute calculations of the 
spectra were made using version in ENDF/B cross sections and an S n discrete ordi- 
nates multigroup transport code. 

Generally excellent agreement was obtained for Be, CHg, Pb, and Mo, and good 
agreement was observed for Nb, although discrepancies were observed for some energy 
ranges. Poor comparative results, obtained for Ta and W, are attributed to unsatis- 
factory nonelastic cross sections. The experimental sphere leakage flux spectra aren 
tabulated and serve as possible benchmarks for these elements against which reevalu- 
ated cross sections may be tested. 


INTRODUCTION 

An effort continues to develop improved neutron cross section libraries (ref. 1). 

The absolute precision required by reactor and shield analysts is dependent on the sensi- 
tivity of performance criteria to uncertainties in cross sections. Because of the spectral 
and importance variations in reactors, these sensitivities are difficult to specify in gen- 
eral terms and existing problems are recognized (ref. 2). Except for isolated isotopes 
such as hydrogen (H) or carbon (C), which qualify as absolute standards for total cross 
section, and a few isotopes such as lithium-6 ( Li), boron- 10 ( A B), gold-197 ( A Au), 



238 

and uranium-238 ( U), which ultimately may qualify as secondary standards for cap- 

ture cross sections, an atmosphere of controversy and uncertainty still exists (refs. 3 
to 6). 

Cross sections for most of the important elements are gradually being remeasured 

and reevaluated. For example, the state of the art and the dilemmas facing the cross 
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section evaluator are well documented by Davey (ref. 7) for the case of U fast neu- 
tron data. In addition, a discussion by Orphan (ref. 8) has reviewed cross sections for 
reactors and shields and has drawn attention to recent developments in techniques for 
measuring elastic and inelastic scattering cross sections. For example, improved fast 
neutron time -of- flight LINAC based systems have permitted measurements for many 
elements up to ~ 9 MeV (ref. 9). Experimental resolution is sufficient to determine in- 
elastic scattering peaks for lower energy isolated levels for middle weight elements and 
to indicate considerable structure for the unresolved higher energy levels. These 
measurements raise questions about the adequacy of representing inelastic cross sec- 
tion as smooth average curves in ENDF/B files. In addition, for some elements, in- 
elastic down scattering is represented crudely by the evaporation model with spectra 
defined by an energy dependent nuclear temperature. 

The situation is perhaps analogous to that for capture cross sections in the kilo- 
electron-volt region where the resolution of the best time- of- flight systems with linear 
accelerators is inadequate to define isolated level structure above a few keV. A knowl- 
edge of the structure is required for fast spectrum reactors and shields. The pitfall of 
assuming that the capture cross sections are smooth is illustrated by recent calculations 
of sphere transmission data (refs. 10 to 13) which account for resonance structure of 
the unresolved cross sections. These Monte Carlo calculations used resonance param- 
eters selected from merged and separated spin statistics and have shown relatively 
large errors in the absolute cross sections at 24 keV obtained from a smooth cross 
section approximation. 

- For neutron flux calculations, the use of modern transport methods and references 
to programs coded for digital computers using the discrete ordinates and the Monte Carlo 
methods have been summarized by Maienschein et al. (ref. 14). What are needed by re- 
actor and shield analysts are nuclear data files of sufficient precision such that designers 
have a high degree of confidence in applying the increasingly available neutron and 
gamma- ray multigroup transport computer codes. The credibility of the design methods 
using these codes that employ ENDF cross section libraries lies in the degree of corre- 
lation obtained between results of integral experiments and their analyses. Integral ex- 
periments for both critical assemblies (ref. 15) and so-called benchmark shield configu- 
rations have been performed (ref. 16). Benchmark experiments consist of well docu- 
mented integral experiments with various media and ideally should be simple in geometry 
and homogeneous in composition. They should be designed to permit measurement of 
absolute neutron energy spectra and angular fluxes so as to facilitate the checking of 
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transport calculational methods and the completeness and precision of neutron cross 
section sets used. A recent review by Straker (ref. 17) discusses the integral experi- 
ments that have been performed and presents a background on the evolution of benchmark 
problems. 

It has been argued that accurate prediction of measured integral quantities such as 
critical masses, material reactivities, and spectral detector indices have inherent 
limitations and may contain compensating effects that mask potentially large uncertain- 
ties in reactor spectra and/or cross sections. An example of the limitations of meas- 
uring spectral indices in spherical samples within a fast reactor assembly as commen- 
taries on nonelastic scattering cross sections has been presented by Davey and 
Amundson (ref. 18). Another complication is the presence of several elements in addi- 
tion to fissionable materials in critical assemblies that makes it difficult to attribute 
discrepancies to specific isotopes. Therefore, alternative integral experiments have 
been sought that have elemental compositions. 

It is clear that experiments that are projected as definitive tests of the precision of 
ENDF/B cross sections should provide differential data so that comparison with calcu- 
lations can serve as commentary on the validity of cross sections over particular energy 
ranges. In this connection the organic liquid scintillator measurements of spectra for 
uncollided fission neutrons transmitted through thick samples of material (refs. 19 and 
20) was a pioneering step in obtaining differential data to test total cross sections. 

Here, comparison of measured and calculated transmission spectra from ~0. 8 to 
12 MeV for many materials served primarily as tests of minima in neutron total cross 
sections. 

Another method for obtaining differential data from idealized experiments employs 
time-of-flight techniques to measure fast neutron angular flux spectra at different posi- 
tions in simple homogeneous assemblies. Representative of these experiments is the 
work of Profio, Cerbone, and Huffman (ref. 21) with relatively thick spheres of paraffin 
and lithium hydride and centrally located photofission sources. Since the sphere is a 
one -dimensional medium and, if homogeneous and composed of one element, it is the 
ideal integral experiment. Neutron energy spectra and angular fluxes are measured by 
extracting a collimated beam from various positions and at various angles in the media. 
Interpretation of the measured angular flux spectra by multigroup transport calculations 
requires (1) a careful absolute experimental method for normalizing the calculated and 
experimental data and (2) a higher order angular quadrature than ordinarily employed to 
calculate scalar flux spectra. The experience of Profio, Cerbone, and Huffman (ref. 21) 
in calculating the measured anisotropic distribution of angular flux, particularly in the 
forward direction, was that it would be necessary to use a very fine angular mesh that 
would require considerable computer storage. The alternative is to use asymmetrical 
quadrature sets to calculate highly forward directed neutron angular fluxes as proposed 
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by Cerbone and Lathrop (ref. 22). However, the angular quadrature used in many reac- 
tor and shield designs is a compromise between angular resolution and computer storage 
requirements. 

The need for clean, homogeneous integral experiments against which multigroup 
transport calculations may be tested is apparent. In the present study, an americium- 
beryllium (Am- Be) neutron source was employed at the center of various spherical 
media. Sphere leakage spectra were measured using an NE213 organic liquid scintil- 
lator neutron spectrometer in the energy range 0. 5 to 12 MeV. Materials studied were 
beryllium (Be), polyethylene (C^), lead (Pb), niobium (Nb), molybdenum (Mo), tanta- 
lum (Ta), and tungsten (W). Preliminary portions of this work were reported earlier 
(refs. 23 and 24). 

It is suggested that these sphere experiments can serve as benchmark data for the 

measured energy range. Future work should be encouraged using other spheres and 
252 

californium -2 52 ( Cf) fission sources and augmenting the energy range of measure- 
ment by employing proportional counters for the keV energy range. 


EXPERIMENT 

A radioactive Am-Be (a,n) source was used in conjunction with a fast neutron spec- 
trometer to measure neutron spectra in the energy range from 0. 5 to 12 MeV. The 
(o!,n) source was placed at the center of spheres of various test materials, and the 
leakage flux spectra measured at a distance of 200 centimeters from the source to the 
spectrometer. The experimental arrangement is shown in figure 1 where both the 
source and spectrometer are located 160 centimeters above the floor. Background con- 
tributions due to neutrons scattered into the spectrometer from the walls of the room 
and from the supporting structures were measured using a paraffin shadow cone. Back- 
ground corrections to the measurements were negligible at the higher neutron energies 
but increased to 30 percent at 500 keV. 


Test Material Spheres 

The spheres used are fairly thick, being generally of the order of several mean 
free paths at 4 MeV so that a large percentage of the source neutrons have interacted 
with the sphere before arriving at the spectrometer. Materials studied were Be, CHg, 
Pb, Nb, Mo, Ta, and W. For the heavier materials little energy transfer occurs by 
elastic processes and for Be the (n, 2n) processes are important. The experiments are, 
therefore, primarily sensitive to the nonelastic cross sections and in particular to the 
spectra of the inelastically scattered neutrons. Since a large percentage of the neutrons 
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from a Be (or, n) source lie in the energy range 2 to 8 MeV, the experiment is particu- 
larly sensitive to inelastic reactions occurring in this energy range. 

Each test sphere employed contained a central cavity 6. 15 centimeters in diameter 
to accommodate the Am- Be source. The metallic spheres used were machined or cast 
of high purity metal with the exception of Nb that contains 0. 95 percent by weight of Zr. 
The spheres containing Mo and W were fabricated of 0. 050- centimeter steel sheet and 
filled with Mo powder or fine W balls. In order to study a thick W case, a second, 
larger W sphere was composed of fine W balls in a 0. 48-centimeter steel container. 

The composition and sizes of the spheres are given in table I. For the case of Pb the 
isotopic composition can vary, depending on the geographical location of the source of 

the Pb. The Pb used in the sphere has a measured isotopic composition close to that of 

208 207 

natural Pb; the measured abundances are 51.37 percent Pb, 20.50 percent Pb, 
26.83 percent and 1.29 percent ^^Pb. For purposes of inter comparison with a 

material of well known cross sections, a CHg sphere made from high density polyethyl- 
ene was included. 


Neutron Spectrometer 

Neutron flux spectra from the Am-Be source and the various spheres were meas- 
ured using a 5- by 5-centimeter cylindrical NE213 liquid scintillator proton recoil spec- 
trometer. The virgin (a, n) spectrum used as input in the transport calculations was 
obtained from measurements using a smaller scintillator. The large liquid scintillator 
has the advantages of high sensitivity and well calibrated neutron response functions 
(ref. 25) while the smaller one has better resolution. Although the resolution of scintil- 
lation spectrometers is poor compared with that of time -of- flight measurements, it has 
been measured and is comparable with the energy widths of the neutron groups used in 
the multigroup transport calculations. A direct and more convenient comparison of the 
experimental spectra with multigroup transport calculations is possible, provided that 
the spectrometer resolution function is folded into the multigroup analytical results. 

This procedure is discussed in the next section of this report. A curve of the energy 
resolution of the 5- by 5- centimeter NE213 liquid scintillator is presented by Straker 
(ref. 20) and indicates the resolution to be 10 percent at 11 MeV and to become progres- 
sively poorer to 20 percent at 3 MeV and 35 percent at 1 MeV. For a cylindrical scin- 
tillator as large as the 5- by 5- centimeter size, care must be exercised in determining 
the neutron spectra from the measured proton recoil spectra due to additional contribu- 
tions caused by multiple scattering and carbon interactions in the scintillator. These 
effects are reviewed in detail by Shook and Pierce (ref. 26). 

In the present experiments the measured proton recoil spectra were reduced using 
the calibrations of Verbinski, et al. (ref. 25) and the unfolding code, FERDOR (ref. 27) 
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for the case of the detector in a parallel beam of neutrons incident on the curved surface. 
The present sphere measurements are in an essentially parallel beam geometry because 
of the large separation distance between the source and detector so that the response 
functions generated for this mode of operation are directly applicable. To assure that 
small differences in scintillator size or spectrometer resolution do not cause a signi- 
ficantly different response matrix, monoenergetic neutron spectra were measured at 
2. 8 and 14. 7 MeV; detailed comparison indicated that response functions agreed satis- 
factorily with the response functions of reference 25. 

Pulse shape analysis and two-parameter data acquisition were used to discriminate 
against y-ray pulses. A block diagram of the spectrometer is shown in figure 2 along 
with the photomultiplier tube base circuit. A linear pulse is taken from the ninth dynode 
of the photomultiplier tube, amplified and sent to one side of the analyzer. A pulse 
whose amplitude depends strongly on the shape of the linear pulse is taken from the 
12th dynode. Commercially available amplifiers and analyzer are used. The Owens 
type pulse shape circuit utilizing space- charge saturation is built into the dynode chain 
of the RCA 8575 tube. The two diodes shown in the circuit eliminate the unwanted large 
negative part of the pulse obtained from this type of circuit (ref. 28). With 2000 volts 
across the tube dynode chain, it is possible to separate neutron and y-ray interactions 
down to 20-keV electron energy, which corresponds to a proton energy of 200 keV. This 
separation has been achieved in a 10-to-l ratio of gamma- to -neutron field. 

The spectrometer has been calibrated using a sodium-22 ( Na) y-ray source. A 
calibration curve is shown in figure 3. The calibration energy of 1. 12-MeV electron 
energy corresponds to the half height of the Compton edge. Additional information on 
the spectrometer and data reduction methods used here is described by Shook and Pierce 
(ref. 26). They also measured neutron spectra for a small polonium -beryllium (Po-Be) 
neutron source with three organic liquid scintillators, two smaller in size than the 5 - 
by 5-centimeter detector, and compared results with the spectra measured by the cali- 
brated 5- by 5-centimeter NE213 detector. The spectra for the smaller cylindrical 
detectors were obtained by differentiation of the proton recoil pulse height spectra. The 
approximations involved in the differentiation analysis were indicated to have small 
effects on the precision of the neutron spectra measured with the smaller scintillators 
but did introduce significant error for the larger 5- by 5-centimeter detector. The ad- 
vantage of using the smallest cylindrical scintillator, which measured 1.2- by 
1. 3 -centimeter, is that its resolution is considerably greater than that of the 5- by 
5-centimeter detector. For example, the measurements reported by Shook (ref. 26) 
showed the smaller detector to have a resolution (widthat half maximum) of 8. 6 percent 
at a monoenergetic neutron energy of 2. 84 MeV compared with a resolution of 20 percent 
for the 5- by 5-centimeter detector at this energy cited by Straker (ref. 20). Similar 
results were obtained for proton recoil data corresponding to a neutron energy of 
0. 5 MeV. Therefore, the small Po-Be source spectra measurements made by Shook 
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(ref. 26) using the 1.2- by 1. 3-centimeter cylindrical scintillator are believed to be one 
of the better resolved spectra for small («, n) Be source reactions and serves as an es- 
timate of the virgin source spectrum. This is discussed further in the following sec- 
tion, which describes the Am- Be source used in the present study. 


Source 

The present sphere transmission measurements were made with a large spherical 

241 241 

Am-Be neutron source containing 54 Ci of Am. The 16.7 grams of Am was inti- 
mately mixed with 66. 8 grams of Be powder, cold pressed, and contained in a net vol- 
ume of 84.3 cubic centimeters. The compacts are doubly contained in two stainless 
steel shells, 0. 140 and 0. 152 centimeter thick, with an outer diameter of 6. 10 centi- 
meter and an inner diameter of 5.46 centimeter. 

The source was calibrated by the manganese bath activation method in January 1967 

O 

and the emission rate was found to be (1. 30±0. 08)xl0 neutrons per second. The calcu- 
lations reported herein are based on the measured neutron emission rate so that the 
comparison of the calculated and the measured neutron leakage fluxes from the various 
spheres are on absolute comparative bases. 

Because of the large physical size of the Am-Be source used, the directly measured 
spectra for the larger source cannot be used as input data for the sphere calculations. It 
was necessary to represent explicitly the uniform source region geometry and its stainless- 
steel containers in the transport calculations. The multigroup calculations are therefore 
for a fixed source and employ a virgin (a,n)Be spectrum that is uniformly distributed 
throughout the source region. The virgin source spectrum was obtained from the scintil- 
lator spectra measurements of a small standard Po-Be source that was approximately 
1.2 centimeter long and 1.3 centimeters in diameter and contained approximately 2 grams 
of Be (see ref. 26). These spectra are shown in figure 4 in which measurements by the 
small scintillator using differentiation analysis are compared with the spectrum measured 
by the calibrated FERDOR method (ref. 25) for the 5- by 5- centimeter scintillator. Com- 
parison of the spectra indicates good agreement with the integral flux above 0. 5 MeV to 
within 4 percent. However, there are clear indications of the better resolution for the 
smaller detector as discussed earlier; the source components near 10, 8, 7, and 3 MeV 
are better resolved by the small detector. 

In the region from 1 to 2 MeV, a valley is indicated that has consistently shown up 
in other measurements of the Po-Be source spectra (refs. 29 and 30). Recent measure- 
ments at Oak Ridge (ref. 31) using liquid scintillators, showed essentially identical re- 
sults for neutron spectra from small Pu-Be, Po-Be, and Am-Be sources and also indi- 
cated a valley in the 1 to 2 MeV region. Since the poorer resolution of all measurements 
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in this energy region cannot adequately define the neutron yields, a virgin (ci!,n)Be 
spectrum has been selected and is shown as the histogram in figure 4. This virgin 
spectrum follows the small scintillator measurements except for the valley in the 1 to 
2 MeV range. Here, quite clearly, the poor resolution conceals the true minima be- 
cause of the higher surrounding yields. It was necessary to adjust yields with the 1- to 
2-MeV region to conform with the known resolution functions in this energy range. 

The precision of the virgin spectrum is indicated by the error bars in the figure. 
These error bars represent ±1 standard deviation and are estimated from the counting 
statistics, the energy calibration uncertainty, and the estimated uncertainty in the scin- 
tillator response functions for the case of unfolded spectra. 

A check on the validity of the virgin spectrum shown in figure 4 lies in the agree- 
ment obtained between measurements of the leakage flux spectra from the large Am- Be 
(o!,n) source and the transport calculations of this leakage flux spectra using the virgin 
(o!,n) spectrum as uniformly distributed input. This comparison is presented later. 


TRANSPORT CALCULATIONS 


Fast neutron transport in the spherical shells of test materials containing the 
spherical Am-Be neutron source is described in terms of the space and energy angular 
neutron flux that satisfies the Boltzmann transport equation. In order to provide an 
accurate description of this neutron leakage flux spectra for the present experiments, 
the time independent multigroup form Of the transport equation in spherical geometry is 
employed. This equation is given by 



- M 2 )<P g (r, M) + c Jr) <p g (r, fi) = S (r) 


QQ 

+ ^ (2Z + 1)P Z (M) ^ <Tg(r,g' ~ S)<Pgt(r) (1) 

1=0 g’ 

where 

g = 1, 2, . . . , G and -1 < p < 1 

The notation used is conventional; in particular, a is the macroscopic total cross 
section for group g and Sg(r) is the isotropic fixed source term representing the neu- 
trons originating from the Am-Be(a,n) reaction. 
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The quantity <p (r, \i) is the angular flux for group g at position r for neutrons 
& 

moving in the cone of directions p to ju + dp where p is the cosine of the angle 
between the neutron direction of travel 0 and the radius r of a sphere. The angular 
flux is expanded in a series of Legendre polynomials as 


<P 


CO 


g (r, m) = ^ ( 2z + 1 ) < ^g(r)P z (M) 
1=0 


(2) 


The coefficients <pl.(r) in the expansion for the angular flux are given by 
& 

<P l g (r) =A jT% g (r,iu)P z (M)dM 1 = 0,1,2,... (3) 

where P^(p)isan £th order Legendre polynomial of argument p. 

The scalar flux $ (r) is given by equation (3) for l = 0; that is, 

$ g (r) s <Pg(r) = ^g( r ’ **) d/i ( 4 ) 

The net neutron current along a radius of a sphere J g (r) is given by equation (3) for 
1=1; that is, 


V r > * 



<P g (r, M)M dp 


(5) 


The CTgfe* — g) are coefficients in the expansion in a Legendre series of 
a s (g' — g> M l)> which is the probability per unit path length of a neutron scattering from 
group g' to group g through an angle where is the cosine of the angle be- 
tween the initial and final directions of motion of a scattered neutron. This expansion 
for o g(g* - g, M l ) is given by 


fr s (g’ ~g>M L ) = 



2Z + 1 
4tt 


a^(g f -g)P;(M L ) 


(6) 


The Zth order group transfer coefficients includes elastic scattering and, if present, 
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inelastic scattering and the (n, 2n) reaction. The inelastic scattering and (n, 2n) reaction 
are usually included only for the l - 0 term. 


Multigroup Solutions 


Given the basic group cross sections for the materials involved in an experiment, 
the source spectrum, and the arrangement of the materials in spherical shells, a solu- 
tion to equation (1) is obtained with the use of a computer. A solution to equation (1) 
consists of obtaining the discrete values of the angular flux <pg(r. , /* m ) for group g at 
position r. and direction that approximate the angular flux <^g(r, \i). This ap- 
proximate solution to equation (1) is obtained by using a code based on the discrete or- 
dinates S n method that is a descendant of the DTF-IV code (ref. 32). Since the angular 
flux is expected to be anisotropic for this geometric arrangement of neutron source and 
spherical shells of test materials (i.e. , the angular flux is expected to be forwardly 
directed), 16 directions were chosen to represent the angular flux. The 16 direc- 
tions (comprising the Sjg calculation) were obtained from the Gauss- Legendre quadra- 
ture set. The adequacy of this Sjg quadrature set was tested by comparing the com- 
puted leakage spectra for the Ta sphere experiment with the results obtained for an 
quadrature set; excellent agreement was obtained between the two calculations with the 
Sjg calculation differing by less than 1 percent for all energy groups from the more 
accurate Sgg calculation. 

Microscopic cross section data for all of the test materials except W were obtained 
from the ENDF/B version HI data files. These data were processed using a version of 
the SUPERTOG code (ref. 33) for the group split shown in table II excluding the thermal 
group (group 67). This 67 energy group split was chosen so that the energy width of the 
group was much less than the resolution of the liquid scintillator spectrometer and 'made 
the cross section averaging procedures insensitive to the assumed averaging spectra. 
Elastic scattering cross sections were included through the l - 3 term of the scatter- 
ing source term on the right side of equation (1). Any inelastic scattering and (n, 2n) 
reaction cross sections were treated through the 1=0 term; that is, these two types 
of neutron scattering were assumed to be isotropic in the laboratory coordinate system. 

For the W sphere the version of SUPERTOG that was available was not able to 
process all of the necessary data files. Therefore, the W data were obtained from the 
GAM-n code (ref. 34). The basis for the data in the GAM-II data files is reference 35. 
Fifty- one energy groups were used with the first 49 being the first 49 energy groups of 
the GAM-II code that covered the energy range from 14.9 to 0. 11 MeV, with each group 
having a lethargy width of 0. 1. The 50th and 51st groups of this energy group split were 
the same as the 66th and 67th of the 67 energy group split. This 51 energy group split 
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is unsatisfactory at the high energy end as the group width is of the order of the spec- 
trometer resolution. 

The thermal group for all of the materials used in this program were obtained from 
the GATHER-n code (ref. 36). Although the neutron groups below 0. 1 MeV are not 
needed for comparison of the calculated neutron leakage spectra with the measurements 
obtained with the spectrometer, they are needed to complete the problem description. 
This permits an assessment of the neutron balance tables for any code malfunctions. 

When performing calculations with the discrete ordinates S n method, the highest 
order l used in the scattering source term of equation (1) defines a P^ calculation. 
The calculations made herein are thus Pg calculations. A calculation for the Ta 
sphere experiment indicated that a lower order P^ treatment of elastic scattering was 
adequate for computing the neutron leakage spectra. A comparison with the neutron 
leakage spectra computed for a Pg treatment of elastic scattering indicated that the 
results for a Pj calculation differed from those of a Pg calculation by less than 
2 percent for all groups. For all of the metallic test materials, similar results were 
obtained. The adequacy of the Pj calculation is due to the dominance of the inelastic 
scattering and, in the case of Be, the (n, 2n) reaction rather than elastic scattering is 
responsible for the removal of neutrons from one energy to a lower energy region. On 
the other hand, a Pg treatment of elastic scattering was necessary for the CHg sphere. 
Therefore, since the scattering order does not appreciably affect computer time, all 
calculations reported herein were performed with a Pg treatment of elastic scattering 
to insure maximum accuracy. 

Smoothed Leakage Spectra 

For direct comparison of the multigroup calculations with the liquid scintillator 
spectrometer measurements of the neutron leakage spectra at 2 meters from the center 
of the test spheres, the multigroup results are smoothed using the known resolution 

Q 

function. The known Am- Be source strength of 1. 30±0. 08x10 neutrons per second is 
used for the calculations. The spectrometer measurements have been reduced by the 
FERDOR code (ref. 27), and the data are given as a differential flux in units of neutrons 
per square centimeter per second per MeV. Two methods of obtaining the calculated 
differential flux at the detector position are feasible. In the first method the flux is 
calculated directly by the S n method from the source to the detector. This method 
results in a large number of spatial mesh points in the intervening region and a corre- 
sponding increase in the computing time. With this method the calculated differential 
flux is given as (method 1) 
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(V 


$(r d ,E) = Yj t?g^ r d^ AE g] P ^ E g - E ^ AE 
g=l 


where 


P (Eg - E) AE 


J~ 2 

> no 


exp 


■(E g - Er 

2a 2 


AE 


( 8 ) 


In equation (7) ^(r^) is the scalar flux for group g at the detector position r d , AEg 

is the energy width of group g in MeV, p(E — E) AE is the probability that an event 
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at energy Eg will yield a response in AE about energy E, and a is the energy reso- 
lution of the spectrometer in MeV corresponding to the full width at half- maximum of 
the Gaussian distribution. In equation (7) the summation extends over G groups where 
the 67 groups may be subdivided into a larger number of groups for this folding proce- 
dure. 

In the second method the neutron differential flux at the detector position is calcu- 
lated from the neutron leakage computed at the surface of the spherical assemblies 
2 

assuming a 1/r spatial dependence; that is, <£(r d ,E) is given by (method 2) 


<Kr d ,E) = 



L g (r s )/AE gj p ( E _ AE 


4nr~ 


g 


(9) 


In equation (9) Lg(r g ) is the number of neutrons in group g escaping from the surface 
of a spherical assembly that has an overall radius of r g . 

For the smaller spheres of the test materials, calculations indicate that the two 
methods for calculating the differential flux at the detector position give essentially the 
same results. For the two larger spheres (Pb and CHg), the two methods differ by 
about 5 percent at energies greater than 5 MeV. 


Sensitivity Calculations 

—An important factor in judging -the utility of a comparison between the present sphere 
experiments and transport calculations is the sensitivity of the comparison to changes in 
cross sections of interest. For the high energy region considered here these are the total 
and nonelastic cross sections and the inelastic down scattering model. Several calcula- 
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tions have been made to estimate these sensitivities. Although they were made for speci- 
fic spheres, the results are expected to correlate well with the thicknesses of spheres in 
mean free paths that are shown in table I. 

For example, a calculation of the effect of increasing the total cross section for the 
Ta sphere by 10 percent resulted in only a 2- to 3-percent reduction in the leakage 
fluxes. On the other hand, a 10 percent increase in the nonelastic cross section while 
keeping the total cross section constant resulted in a 10- to 12-percent reduction in the 
leakage fluxes. Corresponding calculations made for the smaller Mo sphere indicated 
that these relative sensitivities were similar but were reduced approximately in propor- 
tion to the reduced sphere thickness in mean free paths. 

The sensitivity of the comparison to the inelastic down scattering model is more 
difficult to estimate, but the sensitivity of calculated spectra to the evaporation model 
parameter used in the energy dependent temperature for continuum scattering was eval- 
uated for one case. In the evaporation model the energy distribution of inelastically 
scattered neutrons is proportional to 


N(E f - E)« E 


exp 



where E' is the incident neutron energy and 6 is the temperature of the residual 
nucleus approximated by the parameter ^E'/a. 

For the case of the Ta sphere, the parameter a was reduced by 50 percent, and 
the consequent increase in effective temperature resulted in a 3 5- percent increase in 
the leakage flux about the 2-MeV region. This result is illustrated later. 

In summary, then, fixed source calculations were made in spherical geometry 
using 67 energy groups, an S n quadrature order of 16 and a elastic scattering 
order of 3. The calculations are on an absolute basis and are directly comparable with 
the measured differential energy flux at the detector position. The high order of the 
calculation assures that, within the uncertainties in the absolute source strength, the 
density of the various materials, and the dimension of the spherical shells, any differ- 
ences between the calculated and measured spectra can be attributed to inadequacies in 
high energy neutron cross sections. 


DISCUSSION OF RESULTS 

The measured neutron leakage flux spectra at 2 meters from the bare Am- Be source 
and for this source at the center of various spheres of test materials are presented. The 
absolute measurements are compared with the absolute transport calculation of the same 
quantity employing the known source strength and the virgin Be(o!,n) source spectrum 
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uniformly distributed within the large Am- Be source. For the convience of others who 
may wish to compare their calculations with the present sphere measurements, the meas- 
ured results are listed on table m. 

Bare Americium- Beryllium Source 

The comparison of measured and calculated leakage flux spectra for the bare Am-Be 
source in air is presented in figure 5. The measured results are shown by the vertical 
bars that represent ±1 standard deviation that is estimated from the counting statistics, 
the energy calibration uncertainty, and the uncertainty of the FERDOR response func- 
tions. 

In order to facilitate direct comparison of the calculations with the measured neu- 
tron flux per MeV at 2 meters, the calculated multigroup flux per MeV at 2 meters has 
been smeared by the known resolution function of the NE2 13 liquid scintillator and the 
results appear as the solid line. 

The large source used has sufficient Be and Fe to modify the virgin (o!,n) spectrum; 
ENDF/B cross sections have been used for these elements. The effects of Fe are small 
but the contributions of (n, 2n) reactions with Be are calculated to result in a net en- 
hancement of source strength by 4. 2 percent in excess of the distributed virgin source 
input. To satisfy conservation, calculations were therefore based on a distributed vir- 
gin source strength 4.2 percent less than the measured strength of the Am-Be source. , 
The comparison of leakage spectra for the bare Am-Be source is shown in figure 5 and 
indicates excellent absolute agreement with calculations over the entire measured range. 
Quite clearly the resolution function of the scintillator spectrometer has smoothed out 
the prominences of the input spectra but the shapes and magnitudes agree. 


Beryllium 

The most important high-energy reaction in Be that affects downscattering is the 
(n, 2n) reaction. It is the comprehensive work by Perkins (refs. 37 and 38) that presents 
a detailed model for this reaction and provides data for the angular and energy distribu- 
tions for neutrons from the various decay modes of the process. An early integral test 
of these cross sections was presented by Alexander, Fieno, and Ford (ref. 39), where 
the slowing down flux was measured in a water medium that contained an Am-Be source 
-■ enclosed in Be spheres. The superiority of the Perkins data over other sets was demon- 
strated. 

The Perkins models are incorporated in the ENDF/B evaluation for Be (MAT 1154), 
and this has been used in the present calculations. However, in the multigroup cross 
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section preparation by SUPERTOG, the energy-angle correlations available in the eval- 
uated data for the (n, 2n) reaction are neglected. In the present calculations, the effects 
of ignoring and preserving this angular correlation were studied by (1) using ENDF/B 
cross sections as processed by SUPERTOG and (2) calculating and using group cross 
sections and transfer matrices by the AGN-SIGMA code (ref. 40). In the second calcu- 
lation, the energy- angle correlation for the two branches of the Be(n, 2n) reaction via 
the 2.43-MeV level were represented to P^ order rather than Pq order. 

The partial cross sections due to the 2.43-MeV level dominate the (n, 2n) reaction 

above 4 MeV in the Perkins model as shown in figure 6. The majority (91. 7 percent) of 

8 * 

these reactions are assigned as leaving Be in a broad 2. 90- MeV excited state with an 

O 

effective excitation energy in Be of 0. 6 MeV. The smaller fraction (8. 3 percent) goes 

8 9* 

to the ground state of Be. Partial cross sections for the excitation of levels in Be 

at 6. 76 and 9. 1 MeV are also assigned by Perkins. The remaining cross section, which 

is a large fraction of the (n, 2n) cross section at lower energies, proceeds by multibody 

breakup. 

The experimental leakage flux spectrum for the Am- Be source enclosed in the 20- 
centimeter -diameter Be sphere is shown in figure 7 and compared with the two calcu- 
lated spectra. The solid curve is based on the data of ENDF/B (MAT 1154). The data 
agree well above 4 MeV but a persistent underprediction below this energy is noted. 
However, inclusion of the Pj component in the (n, 2n) scattering matrix results in a 
definite improvement between calculation and experiment in this region, and is therefore 
important for fast spectra studies in Be. For this source-sphere configuration, (n, 2n) 
reactions in the beryllium sphere result in an additional calculated net source enhance- 
ment of 26.6 percent. 


Polyethylene 

The results for the polyethylene sphere serve as a partial check of the absolute pre- 
cision of the present experimental technique. The leakage flux spectra at 2 meters are 
shown in figure 8, where measurements are compared with the transport calculations 
using ENDF/B data for H (MAT 1148) and C (MAT 1165). Although the hydrogen cross 
sections serve as a standard and the total cross sections for carbon qualify as standards, 
the nonelastic cross sections for carbon are not as well known. 

The transport calculations for the leakage flux spectra from the surface of this 
sphere to the detector point at 2 meters have been made by the two methods discussed 
earlier. These results are shown in figure 8. Here, the dashed line represents the 
flux directly calculated by the P 3 S 16 discrete ordinate approximation from the sphere 
through the intervening void to the detector position at 2 meters. The solid line repre- 
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sents the leakage flux at the detector position computed from the discrete ordinate 

fluxes at the surface of the sphere by summing over the surface of the sphere and using 
o 

a 1/r variation from the local surface to the detector. It is this method that has been 
used for the smaller spheres measured. 

The two methods of calculation to the detector agree within a few percent over the 
entire range of energy as they may be expected to since the detector point is of the 
order of 10 radii from the surface of the sphere. 

The agreement of measured and calculated absolute flux spectra above 2 MeV is 
excellent. The oscillations in the measured spectra below 4 MeV are not as well follow- 
ed but are averaged quite well by the calculations. The discrepancy that appears below 
1 MeV is significant and may be due to inadequate downscattering matrices for the in- 
elastic partial cross sections of carbon. 


Lead 

The case of the natural Pb sphere serves as a graphic example of the improvement 
in high-energy inelastic cross sections brought about by comprehensive evaluation tech- 
niques. Until the recent evaluation of natural Pb by Fu and Perey (ref. 41) to provide 
ENDF/B data for MAT 1136, the earlier ENDF compilation for Pb designated MAT 43 
represented the only available data to shield calculators and others. The new evalua- 
tion has separated the reasonably well measured total nonelastic cross section into 
analytical partial cross sections. (See figure 9, which is reproduced from ref. 41. ) 

The calculated components and the measured total nonelastic cross section agreed so 
well that the calculated results were adapted by the evaluators. 

The measured and calculated sphere leakage flux spectra at 2 meters from the 
natural Pb sphere is shown in figure 10. A preliminary analysis of these data was pre- 
sented in reference 24. The calculations using MAT 43 are significantly lower than the 
measurements below 7 MeV and rise above the data below 0. 7 MeV. The comparison is 
particularly poor at 2 MeV where the calculated flux using MAT 43 is half of the meas- 
ured flux. On the other hand the calculations using the later evaluation of MAT 1136 
shows considerable improvement below 3 MeV. However, the calculations using 
MAT 1136 result in about a 20-percent underprediction in the energy range from 3 to 
6 MeV, which suggests that further improvement in the recent evaluation for Pb is 
required. 


Niobium 

Nb represents a material of great interest to reactors and for which the cross sec- 
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tions have undergone several revisions since the initial evaluation by Allen and Drake 
(ref. 42). In late 1971 new data suggested the reduction of the inelastic cross sections 
below 1. 5 MeV and above the (n, 2n) threshold and corresponding increase in the (n, 2n) 
cross section. Individual inelastic levels are delineated up to 1.4 MeV with a continuum 
assigned above 1. 4 MeV based on an assigned nuclear temperature. The ENDF/B eval- 
uated cross section set used is updated MAT 1164. 

The measured and calculated leakage flux spectra at 2 meters from the 25- 
centimeter-diameter Nb sphere are shown in figure 11. The calculations using 
MAT 1164 data are in reasonably good agreement with the measurements at the higher 
and lower energies but underestimate the measured fluxes from 3 to 6 MeV and fail to 
represent the apparent structure in the measured spectra around 2 to 3 MeV. 


Molybdenum 

The myriad of isotopes of Mo make evaluation of the average elemental cross sec- 
tions a difficult task. The basis for ENDF/B cross sections for natural Mo is the eval- 
uation of Pennington and Gajniak (ref. 43) designated as MAT 1111. The high-energy 
cross sections of interest here are based on a combination of measurements and optical 
model calculations for elastic scattering angular distributions and on measurements 
of four resolved inelastic levels below 1. 5 MeV and assigned nuclear temperatures for 
the (n, n'), (n, 2n), and (n, 3n) reactions. The evaluation was revised by Pennington in 
January 1972. 

The measured and calculated leakage flux spectra at 2 meters from the Mo sphere 
are shown in figure 12. The agreement between experiment and calculations is gener- 
ally satisfactory and similar in quality to that observed for the bare source. The Mo 
sphere leakage spectra is considerably softer than that of the bare source of course. 

As noted in table I, the Mo sphere is the smallest sphere measured, and the shell thick- 
ness corresponds to but 0. 8 total mean free path at 4 MeV. 

Sensitivity calculations for this shell thickness indicate that a 10-percent change in 
the total inelastic cross section will change the flux 5 percent, which is slightly smaller 
than the experimental error of ±7 percent. On the basis of these present data, we con- 
clude that the ENDF/B data for high-energy cross sections of Mo represents this ele- 
ment very well. 


Tantalum 

Tantalum represents a case for which a new ENDF/B evaluation has become avail- 
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able and for which the present sphere measurements represent an opportunity to com- 
pare two successive ENDF/B sets of cross sections. The present Ta sphere measure- 
ments were first presented in reference 23 and compared with the earlier ENDF/B data 
compiled as MAT 1035 evaluated by Prince (ref. 44) and by Henderson, DeCorrevant, 
and Zwick (ref. 45). 

A comparison of measured and calculated leakage flux spectra at 2 meters for the 
Ta sphere is shown in figure 13. The calculations using the data compilation MAT 1035 
significantly underestimates the measured fluxes for energies less than 4 MeV. Since 
the parameter a used in the evaporation model of the inelastic continuum has a strong 
effect on the spectra, the value was changed from that used in MAT 1035. The modifi- 
cation consisted of changing the parameter a in the evaporation model from 25 MeV” ^ 
to 12 MeV - '*', which resulted in effectively higher temperatures. A significant improve- 
ment is achieved with considerably better agreement below 6 MeV. However, there is 
no assurance that the calculated spectra are improved below 0. 5 MeV, the limit of the 
present measurements. Above 6 MeV the calculations with MAT 1035 data are consis- 
tently higher than the measurements. This is attributed to smoothed inelastic cross 
sections that are lower than the measured isolated inelastic cross sections. 

A recent reevaluation of Ta data (ref. 46), which has been adopted by ENDF/B as 
MAT 1126, has considered the inelastic cross sections in detail. The total nonelastic 
cross sections and the smoothed curves used in the two evaluations are shown in fig- 
ure 14. With the exception of an isolated data point by Owens, all inelastic data points 
lie between the curves. Although MAT 1035 underestimates the data, MAT 1126 appears 
to overestimate it and results in about a 25-percent difference between the two curves 
over most of the high-energy region. 

The third calculated curve in figure 13 was generated using the more recent 
ENDF/B data of MAT 1126. This calculated curve is significantly below the sphere 
spectra measurements at all energies. There are two possible reasons for these dis- 
crepancies. The first is due to the comparison in figure 14 where it is shown that the 
total nonelastic cross sections for the evaluation of MAT 1126 is significantly higher 
than the data. This appears to be the more likely reason for the lower flux spectra 
above 6 MeV calculated for the Ta sphere using MAT 1126. Therefore a smoothed 
curve of the nonelastic data intermediate between the evaluations of MAT 1035 and 
MAT 1126 would satisfy the present sphere spectra measurements. The second reason 
appears to be the poorly chosen value of the evaporation model parameter a. 

These results were communicated to the evaluators of MAT 1126 and a reply by 

Otte witte (unpublished document entitled "Angular Distribution of Neutrons Inelastically 
181 

Scattered from Ta") suggested that a possible cause of the low sphere leakage flux 
spectra that was calculated may have been due to the neglect of anisotropy in inelastic 
scattering. The data of both MAT 1035 and MAT 1126 specify isotropy for inelastic 
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scattering processes. In order to estimate the effect of anisotropy, Ottewitte made 

181 

calculations for the 136-keV inelastic level of Ta, which indicated significantly in- 
creasing anisotropy above ~6 MeV. Unfortunately, these data could not be included in 
the present transport calculations for the sphere. However, calculations for compar- 
able effects of anisotropy have been made for the case of Be as previously discussed. 

In that case it was shown that consideration of anisotropy for partial (n, 2n) cross sec- 
tions did raise the calculated flux spectra but the effects were limited to energies less 
than 3 MeV as shown in figure 7. For the case of Ta this finding would also be true 
since the inelastic reactions in the continuum region result, in general, in large energy 
losses that effectively remove neutrons from the region above ~3 MeV. Therefore, 
consideration of anisotropy for Ta would be expected to result only in an increase in 
calculated fluxes below ~3 MeV and would not greatly affect calculated fluxes above 
this energy. 

It is concluded that, although the newer evaluations (ref. 46) are based on more 
complete and more accurate data than were available for the previous evaluation 
(refs. 44 and 45) and more sophisticated optical model calculations were performed, 
there is still room for improvement. A process of iteration to include the results of 
integral experiments such as those reported herein can lead to more accurate high- 
energy cross sections. 


Tungsten 

Available compilations of cross sections for W and its isotopes are relatively old 
and are in process of reevaluation by members of the Cross Section Evaluation Working 
Group of Sigma Center at the Brookhaven National Laboratory. The evaluation by 
Joanou and Stevens (ref. 35) has been used to interpret the present sphere measure- 
ments. In figure 15 the high-energy data from this evaluation are compared with the 
data evaluated by Prince (ref. 44) and by Henderson, DeCorrevant, and Zwick (ref. 45) 
that presently comprise the ENDF/B files as MAT 1060 (W-182), MAT 1061 (W-183), 
MAT 1062 (W-184), and MAT 1063 (W-186). These isotopic data were updated in 1970 
by Prince. 

Both sets of calculated data employed the individual partial inelastic cross sections 
for the same identified isotopic levels below 1. 5 MeV. Both employed the measured 
elemental total and elastic cross sections for the individual isotopes. Both employed 
optical model calculations at higher energies and an assigned value of the parameter a 
to define neutron emission spectra for the continuum and identical values of a were 
used for all isotopes. Cross sections for (n, 2n) processes were analytically derived 
above the 8-MeV threshold in both cases. The two sets of cross sections for total, 
elastic, inelastic, and (n,2n) cross sections that are compared in figure 15 indicate sig- 
nificant differences in all cross sections. As mentioned before, the verson of 
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SUPERTOG (ref. 33) that was available was unable to process the necessary ENDF/B 
data listings for tungsten. Therefore, the data used were obtained from the GAM n set 
(ref. 34). It is unfortunate that a comparison with the Prince and Henderson evaluation 
was not possible. 

The measured and calculated leakage flux spectra at 2 meters from the two W 
spheres are shown in figures 16 and 17. In both cases a significant and consistent un- 
derestimate of the measured spectra is calculated below 6 MeV; above this energy the 
data are in better agreement. Although the fewer group calculations do not represent 
adequately the source shape at 8 and 10 MeV, this does not affect the spectra at lower 
energies very much. It is suggested that new evaluations of W be compared with the 
present experimental data, tabulations of which are listed in table m. 


SUMMARY OF RESULTS 

A summary and comments on comparison of the measured and calculated sphere 
leakage spectra are presented in the following table: 


Material 

Figure 

ENDF/B MAT 

Comments 

Beryllium 

7 


1154 

Data agree well above 4 MeV. Consideration of anisotropy in 
(n, 2n) reaction component improves results below 4 MeV. 

Polyethylene 

8 


1148(H) 

Excellent agreement above 2 MeV. Discrepancy below 1 MeV 




1165(C) 

attributed to carbon down scattering model. 

Lead 

10 


1136 

Generally good agreement but some under prediction from 3 to 
6 MeV suggests further improvement in down-scattering model 
required. 

Niobium 

11 


1164 

Good agreement above 7 MeV and below 2 MeV. Calculations 
underestimate data at intermediate energies significantly. 

Molybdenum 

12 


1111 

Good agreement over entire range of energy. 

Tantalum 

13 


1126 

Poor agreement over entire range of energy. Discrepancy above 
5 MeV attributed to incorrect inelastic cross section magnitude; 
below 5 MeV to inadequate treatment of continuum and neglect of 
anisotropy of inelastic scattering. 

Tungsten- 182 

16 and 17 

1060 

Early evaluations have not incorporated more recent data and 

Tungsten- 183 



1061 

are in process of reevaluation. GAM II data inadequate below 

Tungsten- 184 



1062 

6 MeV. (ENDF/B data not explicitly used. ) 

Tungsten- 186 



1063 



Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, August 8, 1973, 

503-10. 
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TABLE I. - COMPOSITION AND SIZE OF TEST SPHERES a 


Material 

Form 

Density, 

g/cm 3 

Outside 

diameter, 

cm 

Container^, 

cm 

Shell 

thick- 

ness, 

cm 

Shell thickness, 
total mean 
free paths 
at 4 MeV 

CH 2 

Polyethylene 

0.957 

45.72 


19.81 

4.5 

Be 

Solid metal 

1. 84 

20.32 


7.00 

1.7 

Nb 

Solid metal 

8.4 

25.4 


9.61 

2.0 

Mo 

Powder 

3.73 

22.96 

(0. 05) 

8.38 

.8 

Ta 

Solid metal 

16,60 

24. 13 


9. 02 

3.4 

W 

Fine balls 

11.87 

22. 96 

(0. 05) 

8. 38 

2.0 

W 

Fine balls 

12.26 

44.24 

(0. 48) 

18.29 

4.5 

Pb 

Solid metal 

11.35 

38.0 


15.9 

4.0 


a All materials are chemically pure with impurities in the ppm range. Only the 


Nb sphere has a significant component of 0.95 percent by weight of Zr. 
^Container material, stainless steel. 
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TABLE II. - ENERGY BOUNDARIES OF 67 GROUP SPLIT 


Group 

High-energy 
boundary of 
group, 
MeV 

Energy width 
of group, 
MeV ' 

Virgin source 
fraction 

Group 

High energy 
boundary of 
group, 
MeV 

Energy width 
of group, 
MeV 

Virgin source 
fraction 

1 

11.331 

0.279 

0. 000330 

37 

3. 867 

0. 188 

0. 029836 

2 

11. 052 

.273 

. 001699 

38 

3.679 

. 180 

. 033700 

3 

10.779 

.266 

.004010 

39 

3. 499 

. 170 

. 032400 

4 

10. 513 

,260 

.006221 

40 

3.329 

. 163 

.026700 

5 

10.253 

.253 

. 007526 

41 

3. 166 

. 154 

. 023188 

6 

10. 000 

.247 

. 008642 

42 

3.012 

. 147 

.021269 

7 

9.753 

.241 

. 007721 

43 

2. 865 

. 140 

.018104 

8 

9.512 

.235 

. 005624 

44 

2.725 

. 133 

.015515 

9 

9.277 

.229 

. 004591 

45 

2.592 

. 126 

. 013445 

10 

9. 048 

.223 

. 006684 

46 

2.466 

. 120 

.012154 

11 

8. 825 

.218 

. 010140 

47 

2.346 

. 115 

.009900 

12 

8.607 

.212 

. 012271 

48 

2.231 

. 108 

. 008040 

13 

8.395 

.208 

. 013867 

49 

2. 123 

. 104 

. 007020 

14 

8. 187 

.202 

. 014402 

50 

2.019 

. 192 

. 011401 

15 

7.985 

. 197 

. 016120 

51 

1.827 

. 174 

. 006968 

16 

7. 788 

.192 

. 016170 

52 

1.653 

. 157 

.006287 

17 1 

7.596 

. 188 

. 014484 

53 

1.496 

. 143 

. 006321 

18 

7.408 

.183 

. 012698 

54 

1.353 

. 128 

. 006091 

19 

7.225 

. 178 

. 011839 

55 

1.225 

. 117 

.007730 

20 

7.047 

.174 

. 014082 

56 

1. 108 

. 105 

. 008440 

21 

6.873 

. 170 

. 014806 

57 

1.003 

. 182 

.017775 

22 

6.703 

. 165 

. 014282 

58 

.821 

. 149 

.016573 

23 

6.538 

. 162 

. 013805 

59 

.672 

. 122 

. 015032 

24 

6.376 

. 157 

. 013378 

60 

. 550 

. 100 

. 012457 

25 

6.219 

. 154 

. 013948 

61 

. 450 

.081 

.010163 

26 

6.065 

. 149 

.013832 

62 

.369 

.067 

. 008391 

27 

5.916 

. 146 

.015370 

63 

.302 

.078 

. 009769 

28 

5.770 

.143 

.017200 

64 

.224 

.058 

.007264 

29 

5.627 

. 139 

i .017426 

65 

. 166 

.055 

. 006889 

30 

5. 488 

.267 

.037562 

66 

. Ill 

. Ill 

. 014036 

31 

5.221 

.255 

.038803 

67 

a . 41 4x10” 6 

. 414x10 ” 6 

0 

32 

4.966 

.242 

.037056 




S = 1. 000000 

33 

4.724 

.231 

.038848 





34 

4.493 

.219 

.036369 





35 

4.274 

.208 

.034159 





36 

4. 066 

. 199 

.031177 






a Lower energy limit is 0. 0001 eV. 
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TABLE ID. - MEASURED ABSOLUTE NEUTRON LEAKAGE SPECTRA AT 2 METERS 


FROM SPHERES CONTAINING Am- Be SOURCE 


Energy, 

MeV 

Bare 

source 

Beryllium 

Polyethylene 

(CH 2 ) 

Lead 

Niobium 

Molybdenum 

Tantalum 

Tungsten 

Small 

Large 



o 

Absolute neutron leakage spectra, neutrons/cm -sec- 

MeV 



0.4 



14. 5±5. 5 






280i40 








.5 

44±4 

71±6 

10. 0±4 

160±10 

190±20 

93±8 



165i30 

.6 

40±3 

65±6 

9. 5±2. 5 

150±10 

165±20 

84i6 



140i20 

.7 

39±3 

54±5 

7. 5±2. 7 

135±7 

1 47 ±15 

75i6 



105il5 

.8 

36±3 

48±4 

5. 5±0. 7 

120±7 

110±10 

65i5 

86i9 

92±8 

74il 1 

.9 

35±3 

48±4 

3. 9±0. 7 

115±7 

91 ±9 

61i5 

70±8 

85i8 

70i8 

1.0 

33±3 

48±4 

3 . 9±0. 7 

100±6 

72 ±8 

56i5 

60i7 

75i7 

57i8 

1. 1 

28±2 

47 ±4 

4. 2±0. 7 

95±5 

56±5 

52i5 

52 ± 6 

66i7 

46i8 

1.2 

24±2 

49 ±4 

4. 0±0. 6 

85±5 

48±5 

50 ± 4 

43 ± 3 

58i6 

34i7 

1.3 

24. 5±2 

51±3 

3 . 9±0. 6 

78±5 

39±4 

47 i4 

37 ± 3 

49i6 

29i5 

1. 4 

25±2 

52±3 

3. 85±0. 5 

73 ±5 

29±3 

43i4 

31i2. 5 

41i5 

23. 5 ± 4 

1.5 

23±2 

52±3 

3. 6i0. 5 

65±4 

24±3 ' 

39i4 

24. 5 ± 2. 5 

35i4 

18. 5 ± 4 

1.6 

21. 5±2 

51 ±3 

3. 3±0. 4 

58±3 

24. 5±3 

36. 5 ± 3 

21. 5i2 

29i3 

15. 5 ± 3 

1.7 

21. 5±2 

48±3 

3. 25±0. 4 

55±3 

26. 5±2 

35i3 

19i2 

26i2 

12. 5i2 

1.8 

23±2 

46±3 

3. 55i0. 4 

52 ±3 

26. 5±2 

34 ± 3 

17 ±2 

24. 5i2 

10. Oil. 5 

1.9 

25±2 

44±3 

3. 8±0. 4 

50±3 

24. 5±2 

33i3 

15. 5il . 5 

24. 5i2 

9. OiO. 9 

2.0 

28±2 

42. 5±2. 5 

4. OiO. 4 

49±3 

23 ±2 

32i2. 5 

14±1. 5 

23. 5i2 

8. 3i0. 9 

2. 1 

30±2. 5 

41 ±2. 5 

3. 8±0. 4 

46±3 

21 ±2 

32i2. 5 

12. 5il 

22. Oil 

7. 4i0. 8 

2.2 

3 2 ±3 

39±2. 5 

3. 5±0. 4 

42 ±2 

2 0±1. 5 

31. 5i2. 5 

12. 2il 

19. 5il 

6. Ii0.7 

2.3 

3 2 . 5 ±3 

3 7 ±2. 5 

3. 3±0. 4 

36. 5±2 

17. 7±1. 5 

30. 5 ± 2. 5 

11. 2i0'. 8 

17. 5il 

5. 5i0. 7 

2.4 

33±3 

35±2 

3. 3i0. 4 

33. 0±1. 5 

16. 5±1. 5 

29. 5i2. 5 

10. 5i0. 8 

16. Oil 

5. 3i0. 7 

2.5 

34±2. 5 

31±2 

3. 5±0. 4 

28. 5±1. 5 

15. 3±1. 5 

29. 0i2 

9. 5i0. 6 

15. 7i0. 8 

5. 4±0. 6 

2.6 

35±2. 5 

27. 5±1. 5 

3. 7±0. 4 

25. 5±1. 5 

14. Oil. 5 

29. 5i2 

9. 5i0. 6 

15. 2i0. 8 

5.4i0. 6 

2.7 

35±2. 5 

23±2 | 

3. 8±0. 4 | 

24. 0±1. 5 

15. Oil. 5| 

29. 0i2 | 

9. 5i0. 6 

16. 0i0. 8 

5. Ii0.5 

2.8 

36±2 

21. 5±1. 5 

3. 8±0. 4 

22. 0±1. 5 

15. 2il. 5 

28. 0i2 

9. 5i0. 6 

16. 2i0. 8 

4. 9i0. 5 

2.9 

36±2 

20±2 

3. 65±0. 4 

21. 5±1. 5 

16. 5 ± 1. 5 

28. 0i2 

9. 5i0. 6 

16. 4i0. 8 

4. 9i0. 5 

3.0 

36. 5±2 

18. 5±1. 5 

3. 45±0. 35 

20. 5il. 5 

17. Oil. 5 

2 7 . 5 ± 2 

9. 5i0. 6 i 

16. 4i0. 8 

4. 8i0. 5 

3.2 

3 7 ±2 

18. 5±1. 5 

3. 3±0. 35 

18. 5±1. 5 

17. Oil. 5 

26. 5 ± 2 

9. 4i0. 5 

16. 5 ± 0. 7 

4. 5i0. 5 

3.4 

35±2 

18. 5±1. 5 

2. 9±0. 30 

16. 5±1. 5 

15. 2il. 5 

26. 5 ± 2 

9. 0i0. 5 

16. 5i0. 7 

4. 3i0. 4 

3.6 

35±2 

18. 5±1. 5 

2. 45i0. 25 

14. 5±1 

13. 5il. 0 

26. 5i2 

8. 9i0. 4 

16. 4i0. 7 

4.3i0. 4 

3.8 

37i2 

19±1 

2. 65±0. 25 

13. 2±1 

12. 7il. 0 

26. 0i2 

8. 9i0. 4 

16. 3i0. 7 

4. 2±0. 4 

4. 0 

37±2 

19il 

3. 0i0. 25 

11. 7±1 

12. 7 ± 1. 0 

2 5 . 5 ± 2 

8. 8i0. 4 

15. 7i0. 7 

3. 9i0. 4 

4.2 

37±2 

18. 5±1 

3. 15i0. 25 

10. 2±0. 8 

13. Oil. 0 

25. 5i2 

8. 5i0. 4 

15. 2i0. 7 

3. 9 ± 0. 4 

4.4 

37±2 

18. 5±1 

3. 15±0. 25 

9. 8±0. 8 

13. Oil. 0 

24 . 5 ± 1.5 

8. 4±0. 4 

15. 2i0. 7 | 

3 . 8 ± 0.4 

4.6 

36. 5±2 

18±1 

3. 15±0. 25 

9. 3±0. 7 

12. 8±1. 0 

24. Oil. 5 

8. 4i0. 4 

15. 0i0. 7 

3. 8i0. 4 

4.8 

| 35±2 

17 . 5il 

, 3. 15±0. 25 

8. 6±0. 7 

12. 7±1. 0 

23. 5il. 5 

8. 4i0. 4 

15. OiO. 7 

1 3. 75i0. 3 
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TABLE HI. - Concluded. MEASURED ABSOLUTE NEUTRON LEAKAGE SPECTRA AT 2 METERS 


FROM SPHERES CONTAINING Am- Be SOURCE 


Energy, 

MeV 

Bare 

source 

Beryllium 

Polyethylene 

(ch 2 ) 

Lead 

Niobium 

Molybdenum 

Tantalum 

Tungsten 

Small 

Large 

2 

Absolute neutron spectra, neutrons/cm -sec-MeV 

5.0 

5.2 

5.4 

5.6 

5.8 

6.0 

6.2 

6.4 

6.6 

6.8 

7.0 

7.2 

7.4 

7.6 

7.8 

8.0 

8.2 

8.4 

8.6 

8.8 
9.0 
9.2 

9.4 
9.6 
9.8 

10.0 

10.2 

10.4 
10.6 
10. 8 
11. 0 
11.2 

11.4 
11.6 
11.8 
12.0 

34.5*2 
33*2 
31.5*2 
29*2 
27. 5*2 
25*2 
23.5*2 
20.5*2 
19.5*1.5 
18*1.5 
18*1.5 
17.5*1 
17.2*1 
16.5*1 
16. 0*1 
15.5*1 
14.2*1 
13.7*1 
11.5*0.8 
9.2*0. 8 
7.7*0. 7 
6. 8*0.5 
6. 5*0. 4 
6. 4*0. 4 
6. 2*0. 3 
5. 8*0. 3 
5. 2*0. 3 
4. 2*0. 3 
3. 4*0. 3 
2. 5*0. 3 
1.7*0. 2 
1.2*0. 2 

16. 5*1 
15±1 
13.7*1 
12. 6*0. 8 
12. 0±0. 6 
11. 1±0. 6 
10. 5*0. 5 
9. 5*0. 5 
8. 7*0. 5 
8. 5*0. 5 
8. 0*0. 5 
7. 7*0. 5 
7. 4*0. 5 
7. 0*0. 4 
6. 6*0. 4 
6. 6*0. 4 
6. 1*0. 4 
5. 3±0. 3 
4. 6*0. 3 
3. 8±0. 3 
3. 35*0. 3 
2. 90±0. 25 
2. 75*0. 25 
2. 62±0. 24 
2. 35±0. 2 
2. 20±0. 2 
2. 05±0. 15 
1. 65±0. 15 
1. 35±0. 15 
1. 00*0. 25 

3 . 05*0. 2 
2. 80*0.2 
2. 60±0. 2 
2. 50*0. 2 
2. 40*0. 15 
2. 35*0. 15 
2. 25*0. 15 
2. 20±0. 15 
2. 15*0. 15 
2. 05*0. 10 
1. 90±0. 10 
1. 85±0. 10 
1. 65±0. 10 
1. 52±0. 08 
1. 40*0. 08 
1. 30±0. 08 
1. 20*0. 07 
1. 09*0. 07 
. 98*0. 07 
. 90±0. 07 
. 84*0. 07 
. 80±0. 07 
.78±0. 07 
. 73±0. 06 
. 68±0. 06 
. 60*0. 05 
. 50*0. 04 
. 40*0. 04 
. 31*0. 03 
. 25*0. 025 
. 165*0. 02 
. 125*0. 015 

7 . 8*0. 7 
7. 2*0. 7 
6. 7*0. 6 
6. 3*0. 6 
6. 0±0. 5 
5. 5*0. 5 
5. 0*0. 4 
4. 8*0. 4 
4. 6±0. 4 
4. 2±0. 3 
4. 0±0. 3 
3. 9*0. 3 * 
3. 7*0. 25 
3.6*0. 25 
3. 4*0. 25 
3. 0*0. 2 
2. 75*0. 2 
2. 35*0. 15 
2. 15*0. 15 
1. 90*0. 15 
1. 75*0. 15 
1. 56*0. 10 
1. 45*0. 10 
1. 35*0. 10 
1. 25*0. 10 
1. 15±0. 10 
. 92*0. 09 
| . 75*0. 07 
! . 58*0. 07 
. 42±0. 06 
. 32*0. 05 
. 21±0. 04 
. 10*0. 02 

11.5*1.0 
11.0*0. 8 
10.0*0.8 
9. 5±0. 7 
9. 0±0. 6 
8. 5±0. 5 
7. 7±0. 5 
6. 8±0. 5 
6. 4±0. 4 
6. 2±0. 4 
5. 8±0. 4 
5. 7±0. 4 
5. 5±0. 4 
5. 3±0. 4 
5. 0±0. 3 
4. 8±0. 3 
4. 6±0. 3 
4. 0±0. 3 
3. 5±0. 3 
2. 9±0. 3 
2.65±0.2 
2.40*0.2 
2.30*0.2 
2.25*0.2 , 
2. 05*0.2 
1.85*0.15 
1.50*0.10 
1.27*0.10 
.92*0.08 
.73*0.07 
. 55*0.05 
.39*0.05 
.28*0.05 
.21*0.05 
.16*0.05 
. 11*0.05 

22.5*1.5 
21.5*1.5 
21.0*1.5 
19.5*1.5 
19.0*1.5 
17.0*1.5 
16.0*1 
14.0*1 
12.5*1 
12.0*1 
11.7*1 
11.7*1 
11.4*1 
11.1*1 
10.7*0.8 
10.2*0.8 
9. 4*0. 7 
8. 4*0. 6 
7. 3*0. 6 
6. 5*0. 6 
5. 2*0. 4 
4.7*0. 4 
4. 6*0. 3 
4. 4*0. 3 
4. 0*0. 3 
3. 6*0. 3 
3. 2*0. 3 
2. 9*0. 3 
2. 2*0. 3 
1.7*0. 2 
1.6*0.25 

8.5*0. 4 
7.7*0. 4 
7. 4*0. 4 
6. 8*0. 3 
6.4*0. 3 
5. 9*0. 3 
5. 5*0.3 
5.0*0. 3 
4.5*0. 3 
4. 1*0.2 
3. 9*0. 2 
3. 9*0. 2 
3. 5*0. 2 
3. 4*0. 2 
3. 4*0. 2 
3. 6*0.2 
3. 3*0. 2 
2. 9*0. 2 
2.55*0. 15 
2.2*0.15 
2.0*0.15 
1. 80*0. 1 
1.65*0.1 
1.65*0. 1 
1.55*0. 1 
1. 5*0. 1 
1.4*0. 1 
1.3*0. 1 
1.0*0. 1 

14.5*0.7 
13.5*0.6 
12.5*0.6 
11.6*0.6 
11.0*0.6 
10.0*0.6 
9.0*0. 5 
8. 0*0. 5 
7.2*0. 4 
6. 9*0. 4 
6.9*0, 4 
6. 6*0. 3 
6. 5*0. 3 
6. 5*0. 3 
6. 1*0. 3 
5. 8*0. 3 
5. 3*0. 3 
4. 6*0. 3 
3. 9*0. 3 
3.35*0.2 
3. 0*0. 2 
2. 7*0. 2 
2. 5*0. 2 
2. 5*0. 2 
2. 4*0. 2 
2. 2*0. 2 
2.05*0. 15 
1.85*0.15 
1.45*0. 10 
1. 10*0.08 

3.50*0.3 
3.30*0.3 
3. 10*0.3 
2.90*0.3 
2.90*0.3 
2.65*0.2 
2.50*0.2 
2.30*0.2 
2. 15*0. 15 
2.00*0. 15 
1.85*0.15 
1.80*0. 15 
1.75*0. 15 
1.75*0. 15 
1.70*0.15 
1.65*0. 15 
1.60*0. 1 
1.40*0. 1 
1.20*0.1 
.97*0.09 
. 87*0. 08 
.77*0.08 
.73*0.07 
.68*0.07 
.65*0.06 
.63*0.06 
.56*0.06 
.48*0.05 
.39*0.05 
.34*0.05 
.26*0.05 
.18*0.05 
. 11*0.05 
















- 
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RCA 857 5 



' High Voltage Source 

Figure 2. - Spectrometer block diagram and photomultiplier tube base circuit. 
















67 Group transport calculation using ENDF/B MAT 1148 

and MAT 1165 cross sections and 1/R 2 variation from 
surface to detector 

67 Group direct transport calculation from sphere to 

detector using same cross sections 
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Figure 8. - Neutron leakage spectra for americium-beryllium source 
in polyethylene sphere. 


o(n, x) 

a(n, 2n) 

o(n f n') continuum 

a(n, 3n) 

o(n, n') level 









= 

si 

1 

nwmm 

raMWWiBVMM 

II 

i 

| 

| 

n 

s 

IVfl 

s 


mB 

is 

5 

S 

■a 

Si 

MHMKB 

^B 

Si 


f| 

■ 

jiBies 

m 

g 

1 

■ M 

■5 


= 

m 

= 


= 

§j 

= 

IB 

n 

■ 

nn 

■ 

■■ 

a 

B 

■ 


Figure 9. - Evaluated natural lead nonelastic cross 
sections and their various components (from 
ref. 41.) Natural lead, 53% 208 Pb, 23% 207 Pb, 
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Figure 12. - Neutron leakage spectra for americium-beryllium source 
in molybdenum sphere. 
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67 Group transport calculation using ENDF/B: 

MAT 1126 cross sections ; 

67 Group transport calculation using ENDF/ B_ 

MAT 1035 cross sections (evaporation 
parameter, a, 25 Mev' 1 ) 
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Figure 13. - Neutron leakage spectra for americium-beryllium source 
in tantalum sphere. 
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Figure 17. - Neutron leakage spectrum for americium-beryllium source 
in large tungsten sphere. 
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